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A phase diagram of two Mott-Hubbard planes interacting with a short-range Coulomb repulsion 
is presented. Considering the case of equal amount of doping by holes in one layer as electrons in 
the other, a holon-doublon inter-layer exciton formation is shown to be a natural consequence of 
Coulomb attraction. Quasiparticle spectrum is gapped and incoherent below a critical doping 5 C due 
to the formation of excitons. A spin liquid insulator (SLI) phase is thus realized without the lattice 
frustration. The critical value S c sensitively depends on the inter-layer interaction strength. In the 
t J model description of each layer with the d-wave pairing, 8 C marks the crossover between SLI 
and d-wave superconductor. The SLI phase, despite being non-superconducting and charge-gapped, 
still shows electromagnetic response similar to that of a superfluid due to the exciton transport. 
Including antiferromagnetic order in the t J model introduces magnetically ordered phases at low 
doping and pushes the spin liquid phase to a larger inter-layer interaction strength and higher doping 
concentrations. 
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I. INTRODUCTION 

A realization of atomically sharp interface of a band in- 
sulator (SrTiOa) and a Mott insulator (LaTiOa) and the 
observation of a metallic interfacial layer between them 1 
has prompted a lot of theoretical work on the interfacial 
phenomena when two insulating materials of completely 
different physical origin are joined together 2,3 . In light 
of rapid experimental advances in the growth of digitally 
sharp layered materials, it seems likely that the successful 
construction of atomically sharp interface between two 
(doped) Mott insulators is not very far. In fact, a success- 
ful construction of alternating LCO:LSCO atomic lay- 
ers was demonstrated recently . A well-known fact that 
some high-T c materials contain bilayers of doped Hub- 
bard planes adds to the practical importance of studying 
the bilayer Mott system theoretically. In these regards, 
a simple yet reliable calculation scheme for understand- 
ing the behavior of heterostructure consisting of Mott- 
insulating materials is called for. We provide one such 
theoretical framework in this paper. A closely related 
work considering a correlated bilayer model can be found 
in Ref. 5. 

The situation we have in mind is that of two lightly- 
doped conducting planes subject to Hubbard-f7 repul- 
sion, coupled to each other through short-range Coulomb 
interaction. Each plane is modeled as the Hubbard 
model, or its large-C/ equivalent as the tJ model. In 
particular we focus on the case of symmetric p-n dop- 
ing, with density S of holes in one layer and an equal 
amount of doubly occupied electrons (doublons) in the 
other. The case S = corresponds to decoupled, half- 
filled Mott insulators in each layer. 

The effective holon-doublon attraction results in the 
exciton binding, which pervades the phase diagram and 
results in a number of unusual properties that have 



no analogs in the single-layer model. For weak inter- 
layer Coulomb interaction, phases found for a single-layer 
model have their counterparts in the bilayer system. Due 
to the exciton binding, however, an incoherent regime 
characterized by a gap in the charge excitation and ill- 
defined quasiparticle peaks dominates the low-doping re- 
gion. This phase is christened the "spin-liquid-insulator" 
(SLI) phase. The in-plane electromagnetic response in 
the SLI regime are shown to be influenced by the exciton 
gap; in particular a superfluid response is found within 
each layer even though the overall phase is not that of 
a superconductor. When each plane is modelled as the 
tJ model with long-range magnetic order, the low-doping 
region is shown to be dominated by antiferromagnetism 
and the SLI phase is found for a large inter-layer Coulomb 
interaction and doping concentration, leading to an in- 
teresting possibility to realize a spin liquid insulator state 
without the lattice frustration. 

Since our attempt is quite new, we present the analy- 
sis of the bilayer model in increasing degree of complexity 
starting with the U — ► oo limit of the Hubbard model in 
section II. Basic strategy and formalism for the analy- 
sis of the bilayer physics is laid out. Hints of new phe- 
nomena unique to the bilayer, as opposed to the same 
model in a single layer, are pointed out. In section III, 
we consider each plane described by the t J model with 
the possibility of having a d-wave pairing of the electrons. 
Finally, antiferromagnetic ordering which was ignored in 
section III is included and its results are compared with 
the non-magnetic t J model in section IV. The inter-plane 
coupling is assumed to be the short-range Coulomb in- 
teraction throughout the paper. Single-electron hopping 
across the layers is not considered. We discuss the over- 
all phase diagram and some spectral as well as electrical 
properties of the phases found with particular emphasis 
on those aspects which have no counterparts in the single 
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FIG. 1: (color online) Schematic plot of the two Mott- 
Hubbard bilayer configuration. The chemical potential is 
chosen to have the same amount of holons in one layer as 
the doublons in the other. In-plane excitations in this model 
for doping S < S c (see text for definition of 8 C ) is depicted. 
Single-particle excitation requires a finite energy A e due to 
the exciton gap. In-plane electromagnetic response is that of 
a superfluid even when the whole system is not a supercon- 
ductor. 



layer. 



II. HUBBARD MODEL FORMULATION 

In this section we outline the basic formulation and 
solve the bilayer model when each plane is characterized 
by the Hubbard model with U taken to infinity. This cor- 
responds to the limit where only the constrained hopping 
of electrons is available. Each plane is assumed identical, 
having the same Hubbard interaction U and the kinetic 
energy scale t. The interaction between the planes is 
mediated by a Coulomb repulsion (n a i = c^Cai) 

Hi2 = Vj2( n u-l)(n2i-l) 

i 

for each lattice site i. The two planes are designated by 
a = 1,2. The Hubbard model for plane a is 
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The Hamiltonian for the combined system reads 
H = Hi + H 2 + H12 ■ 



(1) 



(2) 



We consider the case of a symmetric chemical potential 
mismatch /Ltj = ji = —fi2 an d /x > 0, which leads to hole 
doping for a = 1 and the doping by an equal amount of 
doublons in the second layer. 

It proves convenient to first carry out the following 
particle-hole transformation of the doublon (a = 2) layer, 
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to re-write 



H 2 -> -tVc+ 
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Due to the particle- hole transformation for a = 2, both 
layers now appear as hole-doped and the inter-layer 
Coulomb interaction becomes attractive. 

Double occupation of a given site is prohibited in the 
large-t7 limit, and it is convenient to introduce the slave- 
boson representation of the electron operator 



Caicr 



= b + f ■ 



1 - b^bai- 



(5) 



The constraint b^b a i 



J2cr faiafaia = 1 is included via 
Lagrange multiplier X a i . The slave-boson representation 
results in constrained hopping model for each layer: 



H * = -^b+b^f+J^ + ^KAi 
ija i 

+ E ^(baAi + E faiJai* ~ 1) (6) 
i o 

while the inter-layer interaction takes on 



#12 = -v^^y^fe. 



(7) 



The resulting Hamiltonian H = H1+H2+H12 is solved 
in the slave-boson mean-field theory using the following 
self-consistent parameters 
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All order parameters are assumed real and uniform, and 
\ a i = A is also assumed constant, uniform, and equal 
for both planes. Non-zero A b signifies the presence of 
excitons formed of doublon-holon pairs across the layers. 
The mean-field Hamiltonian thus obtained is a sum of 
bosonic and fermionic parts which reads, in momentum 
space, 



H f - 



11° — 

n mf — 



where k 



tion e 



tker 



e kfakafak 
aka 

Y,&akbak - VA b J2(b^btk + b^b 2k ), 



ak 



(9) 



/ _ 

k 



- —k. We introduced the shorthand nota- 
X b £fe + A and e b k = x^ e k + A* + A for the 
renormalizcd fcrmion and boson dispersions, where Ck — 
—2t[cosk x + cos k y ] is the bare band dispersion appropri- 
ate for a square lattice. For later convenience we re-define 
(i in such a way that the boson dispersion becomes 
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e b k = x f (e k + D) + Li 

with D = 4t being the half bandwidth of the bare band. 
After carrying out appropriate Bogoliubov rotation the 
boson Hamiltonian is diagonalized with the quasiparticle 

energy E b = ^ (e\) 2 - (VA^. 

There are two equations one can write down for the 



fcrmion and boson occupation numbers respectively, and 
another three for the three order parameters introduced 
earlier in Eq. (8). In these equations we replace the 
momentum sum J2k by an integral J deD(e) and use the 
simplistic form of the density of states D(e) = 1/(2D). 
This approximation allows us to solve the self-consistent 
equations analytically. Without going through the rather 
straightforward derivations in detail, we quote the final 
results of solving the five equations at zero temperature: 



A = 5 x X D 

2 4 X W 
X f = (1 - <5 2 )/2 



(2 X f D + /i) 2 - (VA b ) 2 



X ' } = WD + xfD){S + \ ] ~ V{A " )2 \(^(^D + t i) 2 -(VA b ) 2 + \j - (VA b ) 2 )^ 



A b = -VA b ln 



i + n + ^(TTri 2 -(vAby 

li + Vm 2 - (VA*>) 2 



(10) 



The solution of Eq. (10) gives two phases, at T = 0, as 
the doping <5 is varied. For 5 < 5 C a featureless state is ob- 
tained with a gap in the single-particle spectrum. There 
is no well-defined quasiparticle peak in this regime. This 
phase is dubbed the spin liquid insulator, or SLI. For 
5 > 5 C the gap closes, quasiparticle coherence recovered, 
and one is in the Fermi liquid (FL) regime. The crit- 
ical doping 5 C which marks the incohcrent-to-coherent 
crossover is zero for V — and grows as V is increased. 
In the following we discuss the derivation of the results 
mentioned above. For convenience we set the overall en- 
ergy scale 2 X j f D to unity. 

From the second of Eq. (10) one obtains 



(1 + 25) V - (VA ) ] = ]p- 25(1 + 5)} 



(11) 



It is clear that self-consistency is only fulfilled for [i > 
VA b , whereas [i — VA b marks the onset of Bose conden- 
sation, or BEC. Exact expressions for [i and VA b can be 
found as 



VA = 



1 + 5 + Se 2 / y 
e 2/v _ i 



2^8(1 + 5)e 1 / v 
e 2 / v - 1 ' 



(12) 



Note that fi > VA b according to (12). The two values 
coincide at the critical doping given by (restoring 2 X ^D 
temporarily) 



The exciton amplitude A b grows as the square root of 
the holc(doublon) density for small doping. The fermion 
part of the mean-field Hamiltonian carries no gap, but 
the boson spectrum has a gap equal to 



iin(^) = \J ii 2 (VA») 2 = - 



5 - 5e 2 ' v 



o2/V 



(14) 



The gap vanishes at 5 = 5 C . 

When 5 exceeds 5 C , Bose condensation guarantees \i 
is fixed to VA b . For 5 > 5 C , the residual boson den- 
sity 5 — 5 C , defined as the condensate fraction, goes into 
the condensate portion of the boson state. The spectral 
weight of the coherent quasiparticle in the 5 > 5 C regime 
grows as 5 — 5 C . The critical doping depends on the in- 
teraction strength V in the manner given in Eq. (13). It 
is zero for V — as expected for a single-layer Hubbard 
model. 

The T = Green's function for the electron in the SLI 
phase, 5 < 5 C , can be worked out 



G e {k,iuj) = ^cosh 2 ^ — 



+ ^ sinh 2 6 q — 



e 'k+q + E Q - iUJ 



(15) 



where 
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Momentum-integrated spectral function A e (uj) is second layer has a gap with the same magnitude as the 
obtained as the imaginary part of the quantity first layer for 8 < 8 C , has metallic conduction, etc. 
J2k G e (fc, iuj ^ oj + i5). Evaluating the sum over k as an 
integral with a constant density of states as before, we 
get the spectral function A e (uj): 

III. T-J MODEL FORMULATION 



± cosh 2 9 q 9(-E b q ~u)J_ D ?>{x b e - E b q - ui)de 

E sinh2 e i *( w - E D f S ^" e + E \ - u °) de - 

(16) 

The step function 9{— w — E b ) appearing in the first line 
is non-zero only if u> < —E b . The step function in the 
second line is non-zero only if w > E h q . Therefore, in 
the range where |w| < min(£^) we get A e {uo) — 0. The 
gap in the electronic spectral function matches that in 
the boson spectrum obtained in Eq. (14), establishing 
8 < S c as the charge-gapped phase. 

Despite the presence of the gap in the single-particle 
sector the SLI phase is not an electronic insulator, but 
a metal. It can be demonstrated by referring to the 
Ioffc-Larkin rule 9 , which states that the conductivity 
er e of the fermion-boson composite system is given as 
(ct 6 )^ 1 = (<7 b )~ 1 + (o-fy 1 . As will be shown in the next 
section, the exciton pairing has a similar influence on the 
electromagnetic response of the bosons as the BCS pair- 
ing on fermions, and the boson conductivity a b in the SLI 
phase is infinite. The fermion part is metallic, however, 
and the overall response of the electron system is that of 
a metal according to Ioffc-Larkin rule. The peculiarity 
arises from the bilayer nature of the model. Extracting 
a charge from an individual layer costs an energy gap 
because one has to break the exciton pair. On the other 
hand, in-plane motion of a holon is always accompanied 
by that of a doublon, which together make up an exciton 
that moves without an energy cost. 

For 8 > 8 C condensation of single bosons allows us 
to write the coherent part of the spectral function as 
A e {uo) ~ (8 — 8 c )Af (lo) and A* (u) is constant near u) = 
as in the non-interacting Fermi liquid. Therefore we 
have a metallic state with coherent quasiparticles for 
8 > 8 C . One can say that 8 C marks an incoherent-coherent 
crossover of the quasiparticles, while the electric trans- 
port property remains metallic for all non-zero doping. 

Analysis for the second Hubbard layer is simi- 
lar. Due to the particle-hole transformation (3) the 
electron Green's function (T r C2i(r)c2 J -) in the origi- 
nal model is mapped onto — (— l) 1 ^^ {T T c^-{r)c2j) = 
— {— l)*~- 7 'G2jj(— t) after the transformation. Likewise, 
the electronic Green's function at (k,w) is obtained 
as G e (Q — k,—u), Q = (it, it), in the transformed 
Hamiltonian. However, the argument in regard to the 
energy gap is done summing over all available momenta 
k, so the distinction between k and Q — k vanishes. The 



Instead of the Hubbard model in Eq. (1) we adopt the 
t J model for each plane: 



C ajcr Caia I 1 ' 



a ^ ^ n ai + J ^ ' S a i ■ S a j . (17) 
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A well-known mapping of the Hubbard model to the t J 
model for large U/t justifies the use of t J model for the 
doped Mott planes. The spin-exchange term introduces 
additional degrees of freedom having to do with d-wave 
pairing of electrons and long-range magnetic order. We 
will deal with the first possibility in this section and con- 
sider the magnetic order in the next. 

Under the particle-hole transformation (3) for the sec- 
ond layer the spin operator (S x , S v , S z ) is transformed 
to (— S x , S y , — S z ), but the inner product Si ■ Sj re- 
mains unaffected. The Hamiltonian, after the particle- 
hole transformation is carried out, is re-written using the 
slave-boson substitution as the mean-field type H m f — 
mf introducing the mean-field parameters 



Hh mf 



X% = (&+•&«>, Xli = EMajJaicr), A b = <&li&2i>, and 
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aij ai 
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The fermion pairing operator is given by P a ij = 

fatffaji - fail fan in H Lf Assuming all mean-field pa- 
rameters are real and uniform and that the fermion gap 
obeys the d-wave symmetry, we can write down 
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H» — 



E € l faker faker ^{[faktf ' nk \ + f^%Jak\\ 
aha ah 
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Al 



Xf(e k + D)+fi 
— 2i(cos k x + cos k y ) 
2JA* (cos k x — cos k y ). 



(19) 
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The self-consistent parameters and x^ follow from The corresponding spectral function contains a factor 



^/ _ j^f (cosfcr - cosfc^) 2 

k ^k 

1 e f 
X* = --^2(cosk x + cos k y )-^j 



(20) 



at zero temperature. The fermion spectrum is diagonal- 
ized with e[ = \J (e{) 2 + (A^) 2 - At zero doping t k is 

reduced to — 2Jx^(cos k x + cosfc y ) and one readily ob- 
tains x f = A^ w 0.33J, as expected from SU(2) sym- 
metry and in agreement with Kotliar and Liu's earlier 
calculation 6 . The hole number in each layer is given by 

Notice that all of the boson-related self-consistent 
equations are identically those of the J = model de- 
rived in the previous section. The only difference arises 
in the numerical value of \^ which enters the overall en- 
ergy scale 2\ j D. For J = we had \ f = (1 — ^ 2 )/2, but 
with J ^ it is determined self-consistently by Eq. (20). 
Following earlier results 6 , we anticipate that x^ wn l re ~ 
main largely constant over a wide doping range. The crit- 
ical density S c is still given by Eq. (13) with appropriate 

■ The incoherent-to-coherent crossover in the electron 
spectral function with the vanishing of the single-particle 
charge gap occurs at this density of holes. The critical 
doping 5 C is plotted in Fig. 2 as empty squares. The 
fermion sector has a <i-wave symmetry gap in the spec- 
trum throughout the whole doping range until A f it- 
self vanishes at a larger doping. Therefore, S c marks the 
separation of the SLI from the nm-dSC (non-magnetic d- 
wave superconductor). The spectral and electromagnetic 
properties of the individual layer is discussed now. 

The zero-temperature electronic Nambu Green's func- 
tion for the region S < 8 C is worked out as 



G e (k,iuj) 



A k+q sinh 2 9 q 



E J 

q ^k+q 
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with Ak and B k given by 



A k = 
B k = 

cos 29k 
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J Af 
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(22) 



5{oj- E f k+q - E\) and 5(w + E } k+q + E b q ), which yields 
zero when ui — 0, provided the boson gap persists. Thus 
the spectral density is zero at ui — for the SLI phase, 
as was the case for J = 0. It should be noted that in 
regard to the single-particle gap, we no longer have the 
pure d-wave symmetry as was the case for a single-layer 
t J model. The existence of exciton pairing renders the 
gap symmetry in the SLI phase essentially that of d + s. 

Now we present the proof that due to the exciton gap, 
the underdoped region 5 < S c exhibits a superfiuid re- 
sponse to an external field. It is clear that the fermion 
sector behaves as a superconductor with a $ = oo because 
of the cZ-wave pairing. It remains to show that the bo- 
son sector also has a b = oo, and we would have <r e = oo 
according to the Ioffc-Larkin composition rule. 

The method for calculating the electromagnetic re- 
sponse to an electric field for a discrete lattice model had 
been developed earlier by Scalapino et al. 8 . Although it 
was used to calculate the response properties of a fermion 
system, the formalism actually applies equally well to 
bosons. We obtain the zero-temperature superfiuid den- 
sity 



Ps(T = 0) = 



(VAH x f ) 2 



E 



sin 2 k x 



sin 2 k v 



(23) 



for the upper or lower plane of the bilayer. It is non-zero 
as long as the exciton pairing amplitude A fc persists. Ex- 
cept for the (A fc ) 2 , the rest of the terms in Eq. (23) 
proved to be nearly independent of doping in our numer- 
ical analysis. As a result p s ~ (A b ) 2 ~ S for small S. 
The exciton pairing is responsible for the superfiuid be- 
havior of the bosons in each layer. The in-plane electron 
response in the non-superconducting SLI phase is that 
of a superfiuid because <j h and a' are both infinite, with 
(T e ~ 5 at low doping. 

We note that non-zero A b = (6ii&2i) implies the total 
electron number in the individual layer is not conserved; 
rather only the sum over the two layers is. The number 
fluctuation within each layer is reminiscent of the similar 
fluctuation in the superconducting ground state and is 
responsible for the similar electromagnetic response. 



IV. T-J MODEL WITH MAGNETIC ORDER 

The previous sections dealt with the physics of strongly 
interacting bilayers in increasing degree of complexity. 
First we treated the case of in-plane hopping subject 
only to the Mott constraint and the inter-layer Coulomb 
interaction that gave rise to the exciton formation. At 
the next level of complexity we introduced the possibility 
of fermion pairing and of superconducting state in each 
layer. In both cases we found an incoherent regime char- 
acterized by a charge gap at low doping that we called 
the spin-liquid-insulator, or SLI. 
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FIG. 2: (color online) Safm (black) and 5™ (red) for differ- 
ent inter-layer Coulomb interaction V/J with fixed t/J = 4 
according to the mean-field decoupling scheme of section IV. 
The empty blue symbols are S c obtained with the magnetic 
moment suppressed to zero (section III). Various phases are 
identified: AFI (antiferromagnetic insulator), SLI (spin liq- 
uid insulator), m-dSC (antiferromagnetic d-wave supercon- 
ductor), nm-dSC (non- magnetic d-wave superconductor). 

In this section, the spectrum of possible phases in 
the bilayer system is enlarged to include the anti- 
ferromagnetic (AFM) order characterized by m, = 
(l/2)(-l) l (/+ T / an - f+J ail ) = m. We consider the ef- 
fects of the magnetic order on the overall phase diagram 
by decoupling the Heisenberg exchange interaction in the 
hopping, pairing, and magnetic channels. In addition to 
S c we have the critical concentration Safm which marks 
the vanishing of magnetic order. There are four phases 
found, depending on the presence/absence of magnetism 
and of the charge gap. SLI phase, which existed for ar- 
bitrary V in the previous analysis, is now only found for 
a sufficiently large V/ J and 5. The low-doping region is 
replaced by the antiferromagnetic insulator (AFI). The 
superconducting phase is also distinguished by the co- 
existence of magnetism as the m-dSC (magnetic d-wave 
superconductor) at low doping and nm-dSC at higher S. 
The four phases collide at a single tetra-critical point. 
Below we present the details of the mean-field calcula- 
tion. 

Using the slave-boson substitution, the fermionic 
mean-field Hamiltonian if£ , becomes 10 

aha ak 

-2Jm5>/+ +QCT / a ^ (24) 

kcj 

where Q is the AFM wave vector (tt, tt). Other notations 
are identical to those used in Eq. (19). The boson mean- 
field Hamiltonian remains unchanged from Sec. III. Di- 
agonalization of the fermion Hamiltonian and working 



out the self-consistent equations for all the parameters 
involved in the Hamiltonian can be done by straightfor- 
ward manipulation 10 . The self-consistent equations ob- 
tained can be solved numerically for each doping S and 
a given inter-layer interaction strength V. As is well 
known from previous mean-field studies 10 , the magnetic 
order vanishes at around 18% doping for a single layer t J 
model with t/ J = 3. In our numerical solution the crit- 
ical doping concentration for the vanishing of magnetic 
order, Safm, does not vary greatly with V. 

The two critical doping concentrations, Safm an d 
S c , are worked out for varying strengths of V in Fig. 2 
with t/J = 4. The value of S c obtained in the presence 
of magnetic order deviates from that obtained in the 
previous section, when magnetic order was suppressed. 
We use 6™ and <5™ m to differentiate the magnetic and 
non-magnetic critical concentrations. Both are plotted in 
Fig. 2. There are two new phases which were unobserved 
in the non-magnetic model, due to the co-existence of 
magnetism. The charge-gapped phase, which was SLI 
in the non-magnetic model, is either an AFI or SLI de- 
pending on whether the magnetic moment is non-zero or 
not. The d-wave superconducting phase also splits into 
non-magnetic (nm-dSC) and magnetic (m-dSC) d-wave 
superconducting phases. The magnetic phase is seen to 
dominate the low-doping region as expected. The spin 
liquid phase is pushed to high-U/ J, high-<5 region due to 
a competition with antiferromagnetic ordering, but still 
has a finite range of existence. The in-plane electromag- 
netic response remains that of a superfluid throughout 
the entire phase diagram, except in the highly doped re- 
gion, not shown in Fig. 2, where A 6 or A^ itself vanishes. 



V. SUMMARY AND OUTLOOK 

In this paper we worked out the phase diagram and 
some of the unique physical properties in a model of 
coupled Mott-Hubbard bilayers. While all the phases 
of a single layer model have their counterparts in the bi- 
layer model for a sufficiently weak inter-layer interaction 
strength, there also exists a new phase with a gap in 
the single-charge excitation and incoherent quasiparticle 
spectra for sufficiently strong V and/or 5. This spin liq- 
uid insulator phase is to be distinguished from either the 
spin liquid obtained by frustrating the magnetic order 11 , 
or the metallic spin liquid state in the high-T c phase di- 
agram. The uniqueness of the SLI phase we find arises 
from the fact that it is a charge insulator even though 
the doping is finite. 

On the other hand, bilayer exciton formation gives rise 
to the superfluid electromagnetic response in each layer 
in the charge-gapped, spin-liquid phase. Such dichotomy 
of the single-particle and two-particle response functions 
can be explored in future experimental setup. 

It will be of importance to establish the existence of 
the exotic SLI phase found in the present paper beyond 
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the mean-field level and applying techniques other than 
slave-particle approaches. Effects of gauge fluctuations 
on the mean-field phases will be considered in the future 
work. Employing dynamical mean-field theory can help 
carve out the phase diagram of the model Hamiltonian 
used in this paper beyond the slave-particle mean-field 
theory. 



ley and Tiago Ribciro and Alex Seidel for collaboration 
on a closely related project 5 . This work was supported 
by Korea Research Foundation through Grant No. KRF- 
2005-070-C00044. 



Acknowledgments 

H. J.H. thanks professor Dung-Hai Lee for initiating the 
discussion of the bilaycr physics during his visit to Berkc- 



* Electronic address: hanjh@skku.edu 
t Electronic address: cljia@skku.edu 

1 A. Ohtomo, D. A. Muller, J. L. Grazul, and H. Y. Hwang, 
Nature (London) 419, 378 (2002). 

2 Satoshi Okamoto and Andrew J. Millis, Nature (London) 
428, 630 (2004); Phys. Rev. B 70, 075101 (2004). 

3 J. K. Freericks, Phys. Rev. B 70, 195342 (2004); Z. S. 
Popovic and S. Satpathy, Phys. Rev. Lett. 94, 176805 
(2005). 

4 I. Bozovic et al. Nature (London) 422, 873 (2003). 

5 Tiago C. Ribeiro, Alexander Seidel, Jung Hoon Han, and 
Dung-Hai Lee, cond-mat/0605284. 



6 Gabriel Kotliar and Jialin Liu, Phys. Rev. B 38, 5142 
(1988). 

7 Menke U. Ubbens and Patrick A. Lee, Phys. Rev. B 46, 
8434 (1992). 

8 Douglas J. Scalapino, Steven R. White, and Shoucheng 
Zhang, Phys. Rev. B 47, 7995 (1993). 

9 L. B. Ioffe and A. I. Larkin, Phys. Rev. B 39, 8988 (1989). 

10 The details of the mean-field decoupling scheme can be 
found, for instance, in Jung Hoon Han, Qiang-Hua Wang, 
and Dung-Hai Lee, Int. J. Mod. Phys. B 15, 1117 (2001). 

11 P. Chandra and B. Doucot, Phys. Rev. B 38, 9335 (1988). 



